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Abstract 

We report a first-principles study on electronic structures of the deformed armchair graphene 
nanoribbons (AGNRs). The variation of the energy gap of AGNRs as a function of uniaxial strain 
displays a zigzag pattern, which indicates that the energy gaps of AGNRs can be effectively tuned. 
The spatial distributions of two occupied and two empty subbands close to the Fermi level are 
swapped under different strains. The tunable width of energy gaps becomes narrower as increasing 
the width of AGNRs. Our simulations with tight binding approximation, including the nearest 
neighbor hopping integrals between tt- orbitals of carbon atoms, reproduce these results by first- 
principles calculations. One simple empirical formula is obtained to describe the scaling behavior 
of the maximal value of energy gap as a function of the width of AGNRs. 

PACS numbers: 73.22.-f, 73.61.Wp 
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Nanoscale carbon materials including fullerenes and carbon nanotubes have attracted a 
great deal of research interest owing to its versatile electronic properties.- 1 ^ Among them, 
graphene fabricated by Novoselov et al. firstly has been studied extensively- Many in- 
teresting properties of this kind layered two-dimensional carbon nanostructure, such as 
the Landau quantization,- the integer quantum-Hall effect and the quantization mini- 
mum conductivity,- have been investigated by several experimental and theoretical research 
groups. Now much attention has focused on graphene nanoribbons (GNRs) with various 
widths, which can be realized by cutting the exfoliated graphene, or by patterning graphene 
epitaxiallyr^ 1 ^ The edge carbon atoms of graphene ribbons have two typical topological 
shapes: namely armchair and zigzag. All zigzag graphene nanoribbons (ZGNR) are metallic 
due to a localized state at the Fermi level. This has been confirmed by scanning tunneling 
spectroscopy and atomic force microscopy— It originates from a gauge field produced 
by lattice deformation. Such a localized state, however, does not appear in AGNRs. Fujita 
et al. have calculated the energy band structure for ZGNRs and AGNRs by using tighting- 
binding approximations (TBA) for the 7r-states of carbon.— They have found that when the 
width (W) of graphene ribbon is 3n-l, where n is an integer, AGNR is metallic; otherwise it 
is semiconducting. However, the first-principles calculations have shown that the hydrogen 
passivated AGNRs always have nonzero and direct band gaps at the local (spin) density 
approximation level.— The energy gaps (E s ) of AGNRs as a function of ribbon width are 
classified into three families, in which E s (W=3n+l)>E 9 (W=3n)>E 9 (W=3n+2).^ 

The capability to control GNRs' electronic properties are highly desired to build future 
nanodevice directly on GNRs. For example, the electronic structures of AGNRs can be 
altered through the chemical edge modification.— Another possible effective way is to ap- 
ply external strain, since previous studies have indicated that the uniaxial strain affected 
significantly the electronic properties of nanoscale carbon material.— Existing theoreti- 
cal works focuse on electronic structure and magnetic properties of GNRs,— ^ nevertheless, 
more attention should be paid to the geometric deformation effect. In this paper, we perform 
ab inito calculation about strain effect on the electronic structure of AGNRs with various 
widths. Theoretical results show that the energy gaps can be tuned effectively by external 
strain. The tunable width of energy gap decreases when the width of AGNRs increases. The 
nearest neighbor hopping integrals between tt- orbitals of carbon atoms are responsible for 
the variation of energy gap under uniaxial strain. 
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We apply density functional theory with generalized gradient approximation (GGA) im- 
plemented with the DMol 3 package.— The Becke exchange gradient correction and the Lee- 
Yang-Parr correlation gradient correction are adopted.— The basis set consists of the double 
numerical atomic orbitals augmented by polarization functions. The calculations are all- 
electron ones with scalar relativistic corrections. Self-consistent field procedure is carried 
out with a convergence criterion of 5.0 xlO -5 atomic units (a.u.) on the energy and electron 
density. Geometry optimizations are conducted with convergence criterions of 5xl0 -3 on 
the gradient, 5xl0~ 3 on the displacement, and 5xl0 -5 a.u. on the energy. Medium grid 
mesh points are employed for the matrix integrations, the real-space global cutoff radius of 
all atoms is set to be 5.5 A, and uniformly 22 K points along the one dimensional Brillouin 
zone are used to calculate electronic structures of the AGNRs. 

Here, AGNRs with widths W=12, 13, and 14 are chosen to represent three typical families 
(corresponding to 3n, 3n+l, and 3n+2, respectively), similar to the previous theoretical 
study.~ As an example, the schematic of an AGNR with width W=13 is shown in Figure ??. 
To avoid the effects of the a electronic states near the Fermi level, the dangling bonds of 
edge carbon atom are saturated by one hydrogen atom. All atomic positions of AGNRs 
atoms are allowed to relax by using a rectangular supercell, in which AGNR is set with its 
edge separated by at least 10 A from neighboring AGNRs. As a benchmark, the electronic 
structures of AGNRs without geometric deformation are calculated. The calculated energy 
gaps are 0.55, 0.90, and 0.19 eV for the AGNR with width W=12, 13, and 14, respectively, 
which reproduce the previous DFT results.— 

To investigate the electronic structures of these uniaxial deformed AGNRs, the defor- 
mation of AGNRs is quantified by the strain (e) defined as e=(r-r )/r , where r and r 
(r =4.287 A) is the deformed and initial equilibrium lattice constant along the axial direc- 
tion of AGNRs, respectively— The electronic structures of the deformed AGNRs with three 
different widths are all calculated. The band structures for AGNR with width W=13 under 
five different uniaxial strains are shown Figure ??(a), where £=-4.0%, -0.8%, 3.0%, 7.3%, 
and 10.0% labeled with A, B, C, D, and E, respectively. Note that B and D correspond 
to two turning points (the maximal and minimal energy gaps), while A, C, E are three 
intergradation points. Clearly, they exhibit direct band gaps at T point for all cases. The 
obvious difference among band structures of the deformed AGNRs under five different given 
e values is the positions of two upmost valence subbands (vl and v2) and two lowest con- 
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duction subbands (cl and c2) relative to the Fermi level. When the applied strain (e) is set 
to be -4.0 %, all subbands (vl, v2, cl, and c2) are separated at T point. At the maximal 
energy gap point, we observe that two valence subbands (vl and v2) and both conduction 
subbands (cl and c2) degenerate when the geometric deformation is about -0.8%. These 
subbands are separated again when further elongating the AGNR up to 3.0 %. When the 
strain increases continually to 7.3 % (corresponding to the minimal energy gap point), the 
subband vl shifts upwards while cl moves down to the Fermi level, leading these two sub- 
bands almost degenerate. When e further increases to 10.0 %, four subbands are separated 
again. 

We further investigate their electronic properties and observe an interesting phenomenon. 
The spatial distributions of these subbands (vl, v2, cl, or c2) of the deformed AGNRs with 
given e are plotted in Fig. 2 (c)-(d). For A (e=-4.0 %) case, the valence subband vl is mainly 
contributed by the parallel axial bonds, while v2 is featured by the vertical axial bonds. At 
the turning point B (e=-0.8 %), it is clear that two subbands (vl and v2) degenerate at 
T point and their spatial distributions are swapped as seen in Fig. ?? (b). As increasing 
the strain to 3.0 %, the spatial distributions of cl and c2 are interchanged as illustrated in 
Fig. ?? (c). Comparing with spatial pattern at D (e=7.3 %), this phenomenon happens to 
two subbands cl and vl at E (e=10.0 %) as shown in Fig. ??(d). 

To display more clearly, the variations of energy gap of AGNRs with width W=12, 13, and 
14 as a function of e are shown in Figure ??(a) with filled square, circle, and triangle symbol 
lines, respectively. The calculated maximal values of E g for the AGNR with widths W=12, 
13, and 14 are 1.07, 1.00, and 0.97 eV appearing at e=5.0%, -0.8%, and 9.5%, respectively, 
while the minimal values of E g are 0.02, 0.03, and 0.03 eV, which occurs at £=-4.5%, 7.3%, 
and 1.3%. Although the exact semiconductor-to-metal transition does not achieve by both 
elongating and compressing the AGNRs, it is clear that the value of energy gap is sensitive 
to the applied strain (e). In other words, the uniaxial strain strongly affects the electronic 
structures of AGNRs. This implies that AGNRs can be used to design as strain sensor. It 
is interesting to note that the shapes of calculated curves display zigzag feature for three 
different ribbon widths. The energy gaps change almost linearly between two neighboring 
turning points by changing the e. 

The variations of the energy gaps of three family structures with different widths (W=3n, 
3n+l, and 3n+2, where n=4, 5, and 6) as a function of e are shown in Fig. ?? (b), (c), and 
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(d), respectively. Obviously, there exist following four main common features for all AGNRs. 
(1) The zigzag feature is observed for the deformed AGNRs with large width; (2) The energy 
gap decreases when the width of AGNRs increases without geometric deformation; (3) The 
minimal energy gap of all deformed AGNRs is several meV, while the maximal energy gap is 
sensitive to the width of the deformed AGNRs and its value reduces as increasing the width 
of AGNRs. For example, for the deformed AGNRs with width W=3n (n=4, 6, and 8), the 
maximal values of energy gaps are 1.07, 0.74, and 0.56 eV, which appear at £=9.5 %, 6.6 %, 
and 4.8 %, respectively, as shown in Fig. ?? (b). (4) Clearly, the distance between between 
two turning points becomes shorter when the width of AGNRs increases, which suggests 
that the tunable window of energy gap becomes narrow for the wider AGNRs. 

In general, only including a constant nearest neighbor hopping integral (t) between ir- 
electrons TBA results of AGNRs are different from these by first-principles calculations.— 
However, it could reproduce the DFT results of hydrogen passivated AGNRs through in- 
troducing of an additional edge hopping parameter.— 1 ^ Here, to capture a clearer picture, 
the electronic structures of the deformed AGNRs are also calculated using the TB model. 
According to the geometric optimized results, we find that there are four kinds of carbon- 
carbon bond length in the deformed AGNRs, as labeled by a n (n=l to 4) in Fig. ??, where 
ai and a 2 stand for the inner C-C distances, while a 3 and a 4 describe the edge C-C separa- 
tions. The variation of four kinds of bond lengths as a function of the strain e are shown 
in Fig. ??(a) for the deformed AGNRs with width W=13. Similar results are obtained 
for AGNRs with width W=12 and 14. It is clear that the C-C separations change almost 
linearly with the strain and the deformation leads to the largest change of the inner C-C 
bond length (ai). The change of C-C distance results in variation of the hopping parameter 
between two neighbor carbon atoms (t) change in the deformed AGNRs. For simplicity, 
we assume that the change of t (At) is proportional to Aa linearly. Comparing with the 
change of ti (Ati), the changes of other three hopping integrals are set to At2=0.15Ati, 
At3=0.40Ati, and At4=0.13Ati, respectively. The coefficient before Ati is determined by 
the relative slope as shown in Fig. ?? (a). Four initial hopping parameters for the equi- 
librium AGNRs without deformation are set to be: ti=-2.7 eV, t 2 =-2.65 eV, t 3 =-3.2 eV, 
and t 4 =-2.75 eV. The changes of energy gaps of AGNRs versus ti are shown in Fig. ?? 
(b). Clearly, the TBA results reproduce the main feature of DFT calculations. This result 
shows that the electronic structures of AGNRs can be described by introducing the hopping 
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parameter ti in TB scheme. The change of the hopping parameters of the deformed AGNRs 
are responsible for the variation of energy gaps. 

Recently, Han et al. have measured the size of the energy gap of graphene nanoribbons 
with various widths from 10 to 100 nm.— They found that the energy gap scales inversely 
with the ribbon width. Here, we extend the TBA calculations to get the variation of maximal 
value of Eg as a function of the width of AGNRs (W), as shown in Fig. ?? (c). It is clear 
that Eg max decreases smoothly as increasing the width of AGNNRs, which is independent of 
the family structures. By fitting the calculated curve, we obtain a simple empirical formula, 
Eg ma x (eV)=14.06/W, as plotted in Fig. ?? (c) with the green line. We find that this scaling 
relationship can be used to calibrate the maximal energy gap of AGNRs with large width. 
For example, the maximal value of energy gap of AGNR with W=100 (about 12 nm) is 
found to be 0.14 eV from TB calculation, remarkably, this empirical relation curve gives the 
value of 0.14 eV as well. This result also agrees well with the reported experimental value 
of the GNR without geometric deformation.— 

In summary, the electronic structures of deformed AGNRs are calculated by using ab 
initio methods and TB methods. The energy gaps of AGNRs are predicted to change with 
zigzag shape as a function of the applied strain. The tunable window of the energy gap 
becomes narrower when the width of ANGRs increases. TBA simulations reproduce these 
results by first-principles calculations. We find that the change of hopping integrals between 
7r-orbitals of carbon atoms are responsible for the variation of the energy gap of deformed 
AGNRs. These findings are helpful to construct and design graphene nanoelectronic devices 
in the near future. 
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Figure 1: (Color online) Schematic of an AGNR with width W=13. Here, the one dimensional unit 
cell distance between two dash-dotted lines is represented by r. The blue atoms denote hydrogen 
atoms passivate the edge carbon atoms (black dots). Four kinds of carbon-carbon bond lengths are 
labeled with ai to &4. 

Figure 2: (Color online) (a) Band structures of AGNRs with the width W=13. Here, the uniaxial 
strain (e) is set to be -4.0, -0.8, 3.0, 7.3, and 10.0 % and labeled with A, B, C, D, and E, respectively, 
(b), (c) and (c) The spatial distribution of these subbands (vl, v2, cl, and c2) near Fermi level of 
the deformed AGNRs with different strains. 

Figure 3: (Color online) Variation of the energy gaps (E g ) for the AGNRs with width W=12, 13, 
and 14 as a function of the strain (e). Variation of E g as a function of e for three family structures 
with different widths, (a) W=3n, (b) W=3n+1, and (c) W=3n+2 (n=4, 6, and 8.) 

Figure 4: (Color online) (a) Four kinds of C-C bond lengths of the deformed AGNRs with width 
W=13 as a function of e. (b) The Variation of the band gap obtained by TBA method for AGNRs 
with W=12, 13, an 14 as a function of the nearest neighbor hopping integral ti. (c) The maximal 
value of energy gap of AGNRs obtained from TBA results versus the ribbon widths with three 
family structures (W=3n, 3n+l, and 3n+2, n is a integer). 
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Fig. 1 of Sun et al. 
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Fig. 4 of Sun et al. 
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